用强耦合RANS方法模拟旋翼悬停流场

在旋转坐标系下,将Spalart-Allmaras(S-A)一方程湍流模型和Reynolds-averaged NavierStokes(RANS)方程耦合成一个新的RANS方程,并发展了基于多块重叠网格的强耦合RANS求解方法,用于直升机旋翼悬停流场的数值模拟.为了提高计算效率,针对多重网格方法在多块重叠网格上实施的困难,提出了一种基于重叠网格的多重网格实施方法.通过对Caradonna-Tung(C-T)和ONERA 7A旋翼悬停算例验证了发展的强耦合RANS方法和基于重叠网格的多重网格实施方法的有效性.研究结果表明:发展的基于重叠网格的多重网格方法有较高的计算效率,3层网格的加速比约为7.7;强耦合RANS法的计算精度明显高于传统的松耦合RANS方法,特别是在与阻力相关性能参数的预测中,强耦合RANS方法的预测结果更加精确.     

高超声速数值激波失稳的网格依赖性

高超声速强激波的稳定捕捉仍极具挑战性。目前工程计算中普遍应用的数值格式在模拟多维强激波时通常会遭遇明显的激波不稳定现象,且数值格式的激波稳定性对计算网格表现出严重的依赖性。基于矩阵稳定性分析法,对比了具有不同耗散性质的数值格式稳定捕捉激波的能力,分析了空间二阶精度格式的激波稳定性及限制器对激波稳定性的影响。在此基础上,重点探究了计算网格对激波稳定性的影响规律,研究了网格长宽比和畸变角度对激波稳定性的影响。结果显示,在激波附近通过增大网格长宽比或改变网格畸变角度可有效改善激波捕捉的稳定性;相比于增大网格长宽比,改变网格畸变角度提升激波捕捉稳定性的效果更加明显。在此基础上,结合数值耗散分析对高马赫数下数值激波失稳现象的网格依赖性机制进行了探讨。     

面向混合网格高精度阻力预测的熵修正方法

Roe通量差分分裂格式具有耗散小、分辨率高等特点,在亚跨超声速流场模拟中得到广泛应用,但空间离散时需要对Roe平均矩阵的特征值进行熵修正,使用熵修正会增大耗散,影响阻力的求解精度。本文针对非结构混合网格中Roe格式熵修正的特点,通过改进传统的Roe格式Harten-Yee熵修正方法,提出了一种可提高非结构混合网格黏性计算精度的Harten-Yee熵修正改进方法。利用改进后的方法,完成了DLR-F4翼身组合体算例的计算和对比分析,改进后的熵修正方法残差收敛特性与原始Harten-Yee熵修正一致,计算精度提高,计算结果和无熵修正时基本一致,说明改进后的熵修正方法既保留了使用熵修正带来的程序鲁棒性等优点,同时把熵修正对阻力预测精度的影响降到了最低,验证了方法的有效性。采用改进后方法对第3届AIAA阻力会议的计算模型DLR-F6翼身组合体进行了详细的模拟,分析了网格收敛性和雷诺数的影响。结果表明：改进的Harten-Yee熵修正,更加适用于非结构混合网格的黏性计算,计算精准度达到国际同类CFD软件水平,进一步验证了改进方法的可靠性。     

并行化非结构重叠网格隐式装配技术

重叠网格技术广泛应用于复杂外形及多体分离问题的数值模拟,但是超大规模重叠网格装配的鲁棒性和效率依然是实际应用中的瓶颈。基于格心型的有限体积格式,发展了一种并行化的非结构重叠网格隐式装配技术。该技术采用节点的壁面距离作为插值边界的判断准则,并通过物理边界推进确定出活跃区域。此外,基于网格分区策略实现了隐式装配技术的并行化,为了能够适应超大规模的计算网格,在每一个网格子区内单独建立用于查询的交替数字树（ADT）数据结构,整体网格点在各个子区内进行查询,并通过信息通讯接口（MPI）传递查询结果并进行逻辑判断。该并行化的非结构重叠网格隐式装配技术具有较高的自动化程度,通过多体装配算例以及机翼外挂物分离算例,证明了该方法的实用性。     

多块多重网格法及其跨声速转子内流并行模拟

为提高压气机内部复杂流动计算的收敛速度,发展了一种多块多重网格上的可压缩流动并行计算流体动力学(CFD)方法,数值模拟了跨声速转子NASA Rotor 35内部流动及气动性能,重点分析了对接分区方式和内界面通信模式对并行性能的影响。结果表明:对接分区方法和多块多重网格内界面数据处理方法能够保证串、并行CFD方法结果基本一致,且与实验数据吻合良好;并行精度基本不受分区数目和通信模式的影响,并行效率随着相对通信负载的增大而降低。所发展的并行计算方法对跨声速转子气动分析和设计具有一定的可靠性,为开发大型流体机械复杂流动的基础并行算法及软件提供了参考价值。     

不平衡电网三相四开关变换器预测功率控制

为了提高三相四开关变换器并网电能质量,提出一种适用于不平衡电网条件的有限控制集模型预测直接功率控制策略。分析三相四开关变换器的工作机制和电压矢量变化关系,并建立其功率预测模型。使用αβ静止坐标系下的电网电压以及其90°延迟信号计算功率补偿值,设计价值函数,选择最优电压矢量对应的开关状态。该控制策略无需传统的正负序分离控制及锁相环技术,易于实现。仿真和试验结果表明,在不平衡电网条件下,所提出的控制策略能够有效降低并网电流畸变,消除功率波动,提高并网电能质量。     

旋翼滑流对倾转旋翼机飞机模式下气动特性的影响

建立了考虑旋翼桨叶形状、桨叶片数、桨距和转速等因素的等效桨盘模型,该模型中桨叶对气流的作用被等效为时间平均的源项添加到控制方程中,并通过两个算例验证了模型的有效性。对V-22倾转旋翼机飞机模式高速巡航状态下的全机三维流场进行了数值模拟,分析了旋翼流场特征,定量给出了滑流对气动力系数的影响量。计算结果表明:旋翼滑流明显改变了机翼表面的压力分布,使全机升力系数增大,最大升力系数增量为4.6%;0°迎角下滑流使阻力系数减小了29%,而4°迎角后,滑流使得阻力系数增大,在14°时增量达到了32%;俯仰力矩系数的变化量为4.6%。进一步研究发现,滑流对平尾下表面压力系数分布产生较大影响,而对垂尾影响则较小。      

二维欧拉方程的非结构网格自适应方法

给出提高二维欧拉定常解质量的非结构网格自适应方法。网格的生成采用Ｅｅｌａｕｎｅｙ方法,流场的算法为节点型有限体积法,能量的计算采用Ｖａｎｌｅｅｒ格式,自适应采用网格的增加和删减。对流场特性的探测采用密度物质导数的绝对值。最后计算了Ｎａｃａ００１２翼型的跨音速绕流,结果证明自适应效果很好。     

三维复杂几何外形的表面三角形网格自动生成

按点－线－面的层次数据结构,选取几种有代表性的曲线曲面,实现了三维复杂外形的表面描述。采用曲面到平面的保形变换,在平面内由阵面推进法自动生成三角形网格,实现了三维复杂外形的表面三角形网格自动生成。